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ABSTRACT 

We analyse a catalogue of simulated clusters within the theoretical framework of the Spher- 
ical Collapse Model (SCM), and demonstrate that the relation between the infall velocity of 
member galaxies and the cluster matter overdensity can be used to estimate the mass profile 
of clusters, even though we do not know the full dynamics of all the member galaxies. In 
fact, we are able to identify a limited subset of member galaxies, the 'fair galaxies', which 
are suitable for this purpose. The fair galaxies are identified within a particular region of the 
galaxy distribution in the redshift (line-of-sight velocity versus sky-plane distance from the 
cluster centre). This 'fair region' is unambiguously defined through statistical and geometri- 
cal assumptions based on the SCM. These results are used to develop a new technique for 
estimating the mass profiles of observed clusters and subsequently their masses. We tested 
our technique on a sample of simulated clusters; the mass profiles estimates are proved to be 
efficient from 1 up to 7 virialization radii, within a typical uncertainty factor of 1.5, for more 
than 90 per cent of the clusters considered. Moreover, as an example, we used our technique 
to estimate the mass profiles and the masses of some observed clusters of the Cluster Infall 
Regions in the Sloan Digital Sky Survey catalogue. The technique is shown to be reliable 
also when it is applied to sparse populated clusters. These characteristics make our technique 
suitable to be used in clusters of large observational catalogues. 

Key words: galaxies: clusters: general - galaxies: kinematics and dynamics - large-scale 
structure of Universe 



1 INTRODUCTION 

The clusters of galaxies are the largest dynamically relaxed struc- 
ture we observe in the Universe. The knowledge of their dynam- 
ics and mass distribution is strategic in modern cosmology, since 
it constrains different cosmological models (White & Fren k|199T) 
|Bahcall & Cen][T993l |Borgani et aT^[l997) . Several methods are 
currently used to estimate the cluster masses and mass profiles, 
based on the dynamical analysis of the member galaxi es (|Kaiser 



1987) |Reg6s & G eller 1989; Diaferi o'XGeller||1997| |Girardi et 
al.||1998[>, on the observations of the X-ray emitting gas iCowie, 



Henriksen, & Mushotzkyl|1987| |Borgani et al.||2004^, an d on the 
gravitational lensing effects (Grossman & Narayan 1989, N arayan| 
& Bartelmann 19961. However, the estimation of cluster masses 
remains difficult, mainly because most of the matter in clusters is 
actually dark matter, which cannot be easily observed. 

Within this framework, the dynamics of the cluster outskirts 
is particularly relevant. In the oustkirts, the matter is not yet in 
equilibrium and is affected by an overall infall motion towards the 
cluster centre. As discussed by Cupani, Mezzetti, & Mard irossian| 
(2008 . hereafter CMM08), the dynamics of this' 'non-equilibrium" 
region turns out to be quite well described by the Spherical Col- 



lapse Model (SCM; 


Gunn & Gott|1972||Silk|1974||Peebles|1976| 


|Gunn|1978||Peebles 


1980 1. CMM08 focused on the turnaround ra- 



dius r t (i.e. the radius where the mean velocity of the infalling mat- 
ter balances the Hubble flow) and studied the dependence of the 
overdensity 5 on the radial coordinate r in the region surrounding 
r t . They demonstrated that the large majority of clusters are com- 
patible with a single mass profiles when their mass M is expressed 
in units of the turnaround mass Aft. 

The results of CMM08 provide a description of the overall 
shape of the cluster mass profile in the non-equilibrium region. In 
principle, a way to reconstruct the individual features of these pro- 
files is provided by the relation between the overdensity and the 
velocity of the infall motion in clusters (Regos & Gelle r|1989) : 

Vint 



Winf 



H r 



n u M ° F(S). 



Here Vi n f is the infall velocity (i.e. the bulk peculiar velocity of 
matter towards the cluster centre, defined as positive when directed 
inwards), Hq is the present-day Hubble parameter, S1m,o is the 
present-day matter density parameter and F is an analytical func- 
tion of the overdensity S (CMM08; |Yahil|1985[fvTllumsen & Dayls] 
|1986[ >. If Vi n { and F are known, one can compute the mass profile 
by inverting equation 
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M(r) = -Trpc^o^M.or- 3 [1 + 8(r)] 
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= ^7rpcr,o^M,or 3 {l + F 1 [f2 M °' 6 aj inf (r)] } , (2) 

where /9 cr is the critical matter density, and F^ 1 is the inverse func- 
tion of F. 

To compute the mass profile via equation |2| one must recon- 
struct the overall infall pattern of clusters using the member galax- 
ies as tracers. Such reconstruction may result to be difficult for two 
reasons: 

(i) the radial velocities of member galaxies are affected by the 
presence of local substructures inside the clusters; therefore, the 
values of t>i n f are distributed along r in a blurred band (see also 
CMM08 for details) and 

(ii) the infall velocities of galaxies can not be easily inferred 
from observations, owing to projection effects. 

Nevertheless, in the present work we wish to show that the infall 
velocity approach can be used to estimate the overdensity and mass 
profiles of galaxy clusters. We focus on the dynamical quantities of 
galaxies which can be directly inferred from observation, namely 
the sky-plane distance from the cluster centre, r sp , and the line-of- 
sight velocity, ui os (all these quantities are the moduli of the cor- 
responding vectors). The technique we put forward is based on the 
identification of a 'fair region (FR)' in theredshift space (r sp , uios). 
The galaxies which lie within this region will be proved to be suit- 
able to reconstruct the total matter distribution in clusters via equa- 
tion (2]l. 

We will discuss our technique using a catalogue of simulated 
clusters, to be able to know the overall kynematical properties of 
the galaxy distribution. The results of this analysis are consistent 
with the formulation of the SCM discussed by Cupani, Mezzetti 
and Mardirossian, in preparation (hereafter CMM2010) and hold 
within a range of values of the present-day matter density param- 
eter (0.2 ^ Am.o ^ 0.4, for a spatially-flat Universe). Our mass 
estimation technique relies only on the observable kinematic prop- 
erties of cluster member galaxies, and it is therefore suitable to be 
applied to observed clusters. As an example, we will discuss the 
results of our technique when applied to a subset of clusters from 
the Cluster Infall Regions in the Sloan Digital Sky Survey (CIRS) 
Catalogue (hereafter RD |Rines & Diaferio|2006) . 

In section|2] we describe the simulated data catalogue we use 
for our analysis. In section [5] we discuss thoroughly the statistics 
of the member galaxy velocities. We first provide an empirical def- 
inition of the FR and then we justify this definition using the SCM. 
In section]?] we detail our mass estimation technique and its relia- 
bility. Finally in section[5] we summarize the recipe for estimating 
cluster masses and mass profiles and we draw the conclusion of the 
present work. 



2 THE CLUSTER CATALOGUE 

We performed our analysis on the same data catalogue used by 
CMM08. These data were produced via a large cosmological hy- 
drodynamical simulation of a A cold dark matter (ACDM) Uni- 
verse run by Borgani et al. ( 2004jl, while the catalogue of objects 
(clusters and galaxies) was extracted by |Biviano et al.| ( |2006| >. We 
refer to these papers for the details. 

The simulation was run assuming a flat ACDM cosmology 
with Q.m + = 1, where Qm is the matter density parameter and 
Qa is the cosmological constant density parameter. The adopted 
values of the main parameters (namely the normalization of the 



Table 1. Cosmological parameters of the simulation, compared with the 
mean values (with la uncertainty) obtained from the WMAP, BAO, and SN 
measurements I Komatsu et al. 2008 see text). 
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Figure 1. Statistics of cluster sizes in the simulated data catalogue. Upper 
panel: frequency distribution of the virialization radii r v ; lower panel: fre- 
quency distribution of the virialization masses M v . n c \ us is the number of 
clusters per frequency bin. 

power spectrum as, the Hubble constant h, the present-day mat- 
ter density parameter Qm,o and the present-day density parameter 
of baryonic matter f2bar,o) are listed in Table[T| In the same table 
we also listed the observed values of these parameters, obtained 
by combining the last measurements of the Wilkinson Microwave 
Anisotropy Probe {WMAP) with the distance measurements from 
the baryonic acoustic oscillations (BAO) and from the Type la su- 
pernovae (SN) ( Komatsu et al. 2008). The two sets of values in 
Table[T]evidence a quite good agreement. 

The catalogue we use (Biviano et al. 2006) contains 1 14 clus- 
ters and 963 1 galaxies. The clusters differ in size, with virialization 
radii r v in the range r v = (0.9±2.2)/i -1 Mpc, virialization masses 
M v in the range (8.0 x 10 13 -r- 1.3 x 10 15 ) /i _1 M , and number of 
member galaxies n ga i between 17 and 403. The detailed frequency 
distributions of r v and M v are shown in Fig.^ 

The analysis of dynamics in clusters was performed using the 
cluster catalogue in two ways: superimposing all clusters into a 
synthetic object, and studying all clusters one by one. The former 
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Figure 2. Galaxy kinematics in clusters: C is the cluster centre, G is an 
infalling galaxy and O is a distant observer. 



approach is needed to increase the statistical significance of the re- 
sults, as pointed out by | Vedel"&"H artwick ( 1998), while the latter 
approach is used to test the accuracy of the results when they are 
applied to single-cluster analysis. Throughout the discussion, we 
identify the clusters with the subscript i and the galaxies with the 
subscript g. 

We took into account the particles [dark matter (DM), gas, and 
stellar particles] of the simulation by Borgani et al. J2004) and the 
member galaxies subsequently identified by |Biviano et al.j (20061. 
The particles were grouped into concentric shells in each cluster 
i, using the value of r v; ; as normalization radius, while the mem- 
ber galaxies were attributed to shells according to their radial dis- 
tance from the cluster centre. The shell subdivision is needed to 
make comparable objects of different scale. The shells were de- 
fined using a logarithmical spacing covering the whole radial ex- 
tent from the virialization core to the far outskirts of clusters. We 
considered in our analysis the dynamically relevant quantities ob- 
tained from the simulation, namely the mass of the particles, m, 
and the six phase-space coordinates of both particles and galax- 
ies (x, y, z, v x ,v y , v z ). The position coordinates x, y and z are re- 
ferred to the cluster centre, while the velocity coordinates v x , v y 
and v z are peculiar velocities corrected for the motion of the clus- 
ter centre. We also used the simulated data as 'mock observations', 
i. e. restricting to only two out of six phase-space coordinates: we 
adopted the x direction as the line of sight, and defined the sky- 
plane distance as r sp = \Jy 2 + z 2 and the line-of-sight proper 
velocity as «i os = v x + HqX, The masses of the shells were com- 
puted by integrating the masses of all the enclosed particles. The 
velocities of the shells are the averages of the velocities of all the 
enclosed galaxies. 



3 ANALYSIS OF THE REDSHIFT SPACE 

Fig. [2] represents a cluster C with an infalling galaxy G as seen 
by a distant observer O; in the same Fig. one can see the angles 
e and a used in the following relations. The total peculiar veloc- 
ity vector v g of the galaxy is split into two ortogonal components, 
the radial velocity vector v r - g and the tangential velocity vector 
Vf, g '- v g = v r -, g +Vf. g . Within the widely accepted assumptions of 
spherical collapse model, the overall infall motion of galaxies can 
be considered as isotropic: 



Angle brackets denote the average over all galaxies located at a 
given radial distance r from the cluster centre. It is worth noting 
that (vf. g ) = does not generally imply (vf, g ) = 0; on the con- 
trary, a nonzero tangential velocity modulus is expected in most of 
cases. In fact, as can be seen from the simulations and from the 
SCM (CMM08; CMM10; see also section|3~T), 



^inf 



( v t;g)- 



(4) 



The existence of random tangential motions does not conflict with 
the presence of an overall infall motion of galaxies, since all non- 
radial contribution are averagely suppressed due to isotropy. Here- 
inafter, we will take into account only the moduli of the velocity 
components. 

The relation between the phase-space coordinates r, v r and 
v t and the observational redshift-space coordinates r sp , ui os is ob- 
tained from Fig. fusing trigonometry. Since iti OS : 9 is a proper ve- 
locity, the radial Hubble flow term Hqt must be taken into account. 
Therefore, 



T'sr>:a — Ta COS £ 



{H Q Tg 



I sin e„ + v t .„ cos a„ cos e„ 



(5) 
(6) 



Equations {5} and |6j form an undetermined system and cannot be 
used to compute r g and v r[g when only r sp;9 and ui OS ; 9 are known. 
To overcome this issue, we focus on the subset of all galaxies hav- 
ing r sp;9 ~ r-g and satisfying the condition (|lti OS ;g|) — «inf (see 
below). We will refer to these galaxies as 'tracking galaxies' (here- 
inafter TGs or the TG subset), since they are the best candidates 
to estimate the overdensity and the mass profile using only the ob- 
served redshift-space coordinates. By definition, the TGs have 

(|«Ios ;g (r- sp )|) 



(r) 



(7) 

H r sp 

Substituting equations |5]( and l|7j into equation l|2j, one can com- 
pute the cluster mass profile using only the observable properties 
of the TGs: 

' Mini [»~s P , (|«i os;9 (r ap )|)J 

4 



M(r) 



-7rp cr , fiM,or sp < 1 + F 



.(|«lo B ;s(»" a p)|)' 



H r Bp 



(8) 



The statistical properties of the TG subset are essential to de- 
termine how the infall velocity approach can be used to estimate 
the cluster mass profiles. TGs are not randomly distributed among 
other galaxies, but tend to be gathered in an identifiable region of 
the redshift space. By definition, the TGs have e g ~ 0, which sub- 
stituted into equation |6]) yields ui OS ; 9 — vt- g cosa 9 . Combining 
this result with equation |4j(, one obtains that (a 9 ) ~ 0. Therefore, 
at any given sky-plane distance, the TGs satisfy the condition 



(9) 



(Ko S;9 (r sp )!) ~ ^ {vf,g(r sp )). 

Equation ([5J suggests that TGs are unlikely to be found in the re- 
gion of the redshift space where ui os ~ 0, since both v r and v t are 
generally non-zero. 

These considerations are supported by the analysis of the sim- 
ulated data of our catalogue. We represented in Fig. [3] the de- 
pendence between D r;s and v t -, g (left-hand column) and between 
|itios; 9 1 and Vf, g (right-hand column) as extracted from our data cat- 
alogue. The galaxies are represented with grey dots, while the TGs 
are represented as black dots. The TGs were operatively defined as 
satisfying the conditions 



,(r))=«taf(r), (w* iff (r)) = 0. 



(3) 



' sp;g 



cose > 1 



(10) 
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Figure 3. Dependence between the velocity components of galaxies. The 
grey dots are all galaxies extracted from the simulation. The black dots are 
the TGs, defined as in equations |l0j and in jl 1| with Vtg = 
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Figure 4. Redshift-space distribution of galaxies. The grey dots are all 
galaxies extracted from the simulation. The black dots are the TGs, defined 
as in equations j!0| and j I 1 1 with t^g = 0.2. 
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(11) 



where the parameter r/^ G = 0.2 is used to mimic a quite small dis- 
crepancy around 0. The TG subset selected with this choice con- 
tains 579 galaxies. 

The redshift-space distribution of galaxies extracted from our 
data catalogue is shown in Fig. [4] Also in this case, the galax- 
ies are represented with grey dots, while the TGs are represented 
as black dots. The abscissa is the sky-plane distance in units of 
the virialization radius r S p,v = r sp /r v . Normalization to r v is re- 
quired to make comparable clusters of different size. The ordinate 
is the line-of-sight velocity ui os;g . The resulting distribution shows 
a well-defined trumpet shape, which is typical of caustic surfaces 
and reaches its minimum amplitude at the turnaround (see e.g. |Os-| 
|trikeretal.|1 998 Reg os & Geller|1989| >. The TGs are concentrated 
in two narrow bands of the redshift space, and are almost absent in 
the region where ui os — 0. 




Figure 5. Distribution of galaxy radial velocities, compared with the pre- 
diction of the SCM (from CMM10). The grey-scale represents the number 
density of galaxies for unit surface. The lines are the theoretical radial ve- 
locity profile computed with the SCM, assuming J!m = 0.3 (thick solid 
line) ans Qm equal to 0.2 and 0.4 (lower and upper dashed line, respec- 
tively). 



3.1 The Fair Region of the redshift space 

We can roughly define as 'FR' the region of the redshift space 
where the TGs are mostly concentrated. In fact, an effective defini- 
ton of the FR is required to estimate the mass profile of clusters 
with the member galaxies via the infall velocity approach. We will 
obtain such definition combining a statistical analysis of the sim- 
ulated data catalogue with the theoretical predictions provided by 
the SCM (CMM10). 

According to the condition of equation 191, the FR is defined 
as the region where the mean peculiar radial velocity is almost 
equal to the mean tangential velocity. The v r and v t distribution of 
our galaxies along r is represented in Figs|5]and[6] respectively. To 
determine the grey-scale, we superimposed an orthogonal grid on 
each plane and counted the galaxies in each cell of the grid, in order 
to produce a two-dimensional histogram. The grey-scale represents 
the number density of galaxies for unit surface T>. Both distribution 
are interpolated with the SCM prediction for S1m,o = 0.3 (thick 
solid line) and for Om.o = 0.2 and fijyi.o = 0-4 (lower and up- 
per dashed line, respectively). Despite the notable dispersion due 
to the variance among clusters (see CMM08), the SCM predicts 
quite well the main profile of both v r and v t . The dependence on 
f^M.o is overshadowed by the data dispersion and is generally neg- 
ligible. To determine the location of the FR, we computed the joint 
probability for galaxies to be located in the same grid cell of Figs 
[5] and [6] Fig. [7] represents the distribution obtained by this proce- 
dure (upper panel), compared to the distribution of TGs from the 
simulated data catalogue. Here the grey-scale represents the prob- 
ability parameter V, obtained by normalizing T> to unity over the 
whole plane. The two panels show a remarkable agreement, indi- 
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Figure 6. Distribution of galaxy tangential velocities, compared with the 
prediction of the SCM (from CMM10). The grey-scale represents the num- 
ber density of galaxies for unit surface. The lines are the theoretical tangen- 
tial velocity profile computed with the SCM, assuming S!m = 0.3 (thick 
solid line) ans SJm equal to 0.2 and 0.4 (lower and upper dashed line, re- 
spectively). 

eating that the condition in equation (|9j is suitable to identify the 
FR. 

The empirical definition of the FR is theoretically confirmed 
by the SCM results (CMM10). In a spherically symmetric scenario, 
the FR is theoretically defined as the region where the parameters 
Ar = 1 — r sp /r and Av = 1 — ui os /v r are minimized, in agree- 
ment with the definition of the TGs. Considering equations {3} and 
l|6j, we can write 

Ar = 1 - cose, (12) 
Aw = 1-| M^L _ i J s j ne _|_ j( v cos Q cos £] 

where K v = v t /v r . Equations \12\ and \13\ describe how the in- 
fall velocity profile predicted by the SCM is biased due to projec- 
tion on to the sky plane. This bias depends both on the variance 
in galaxy dynamics (parametrized by K v and a) and on the galaxy 
displacement with respect to the sky plane (parametrized by e). The 
extension of the FR can be therefore determined by taking the ex- 
tremal values of Ar and Av obtained by varying K v and e and a 
around their mean values. In fact, CMM10 showed that the mean 
value of K v is quite close to unity in the non equilibrium region 
for 0.2 ^ SIm.o ^ 0.4, in agreement with equation In partic- 
ular, if Qm,o = 0.3, (K v (r v )) ~ 0.95 for r v ^ 3. So, in order to 
operatively define the FR, we used the following conditions: 

f^K v ^ 0.95iC FR , K:. fk = 1.5; (14) 

- f/FR < e, a r,F R , r? FR = 0.2. (15) 
The parameter K*. FK mimics the variance of K v around its mean 
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Figure 7. Statistical identification of the FR in the redshift space. The upper 
panel is the expected distribution of TGs in the redshift space, computed as 
the joint probability of the distributions in Figs [5] and [6] The lower panel 
is the distribution of TGs in the redshift space extracted from the simula- 
tion. The grey-scale represents the probability density of galaxies for unit 
surface. 

value 0.95, while the parameter ?7i,;fr, mimics a small discrepancy 
of e, and a around their mean value 0. The boundaries of the FR ob- 
tained through equation \\4\ and l |15[ > are shown as dashed lines in 
Fig.[S] superimposed to the grey-scale distribution of TGs (same as 
in the lower panel of Fig.|7|. The agreement between the theoretical 
prediction and the data distribution of our catalogue is remarkable. 
This result shows that the FR, whose existence was originally in- 
ferred from the distribution of the TGs in the redshift space, is due 
to the projection of the infall velocity profile onto the sky plane. 
The SCM is shown to work for 0.2 ^ Om,o ^ 0.4; this allows to 
identify the FR in different cosmologies. 



4 MASS ESTIMATION 

The definition FR in the redshift space provides a way to select the 
galaxies which are suitable to reconstruct the total matter distribu- 
tion of clusters. The galaxies lying within the FR will be referred 
to as 'Fair galaxies' (hereinafter FGs or the FG subset). In this sec- 
tion, we will prove that the infall velocity approach is suitable to 
estimate the mass profiles of clusters using the FG subset. 
Our technique can be outlined as follows: 

(i) identification of the FG subset via the identification of the FR 
provided in Section [3T[ 
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Figure 8. Theoretical identification of the FR in the redshift space. The 
grey-scale represents the probability density of TGs for unit surface (same 
as in the lower panel of Fig.^J. The dashed lines are the boundaries of the 
FR obtained through equation j!4| and equation j!5| with ^*. FR = 1.5 
and r)p R = 0.2. 



Figure 9. Reliability of the raw mass estimate and of the corrected mass 
estimate using the FG subset. The grey dots in the upper panel correspond 
to the values of M TEmig /M g , while those in the lower panel correspond to 
the values of M corr;9 /M g , both represented as a function of rQ^ v;g . The 
average profile and the la uncertainty region of the distribution along the 
radial coordinate are shown in both panels with thick solid lines and dashed 
lines, respectively. 



(ii) detection of the mass profiles, using the redshift-space coor- 
dinates of the FGs. 

We used the simulated data catalogue to test the reliability of such 
approach, comparing the estimated mass profiles obtained with our 
technique with the actual mass profiles of the simulated clusters. 
The simulated data were transformed into 'mock observations' as 
described in Section [2] All galaxies lying within the dashed lines 
in Fig.[8]were regarded as FGs, and their observable redshift-space 
coordinates r sp and u\ os were used instead of the unknown quan- 
tities r and v r to reconstruct the infall pattern and to estimate the 
mass profile of clusters. 

Before applying the mass estimation technique to the FG sub- 
set, we took into account the effect of projection of the galaxy radial 
positions. In fact, while the TGs approximately lie on the sky plane, 
the FGs are in principle distributed all along the line of sight. There- 
fore, using r sp - g instead of r g for the FG subset would introduce a 
bias in the estimation of the mass along the radial coordinate. To 
minimize this bias, we introduce a 'guess' radial distance, defined 
as follows: 



r G;j 



(16) 



where TZ(r sp ) is the mean ratio between r sp;s and r at a given 
sky-plane distance r sp . The guess radial distance rG ;9 will be used 
hereinafter as a replacement of r sp - g when handling the FG subset. 
If the FG distribution is spherically symmetric, 1Z is obtained as 



TZ(r s 



I 



7" 



(7^ 



dri 



v 7 " 



(17) 



■> 

' lew 



dr lo 



where ri os is the line-of-sight projection of the radial coordinate, 
r cut is a cut radius, and u(r) is the number density of galaxies. 
Using the SCM predictions for v(r) and assuming r cut = 7r v 
(since generally we have no appreciable data beyond this distance), 
we approximated equation {XT} with a linear fitting algorithm as 
~R-(fsp) = 5.5 x 10~ 2 f sp + 0.6. This fit is in agreement with the 
data distribution and is used to correct the radial position of FGs in 
our simulated data catalogue. 

The first raw estimate of the cluster mass profiles of is ob- 
tained by substituting the values of ra-, g and u\ OB . g into equation 

M raw;ff = M inf (rG; 9 , |uioa; g\)- (18) 

The function F^ 1 was defined, according to CMM08, as the in- 
verse of the Meiksin approximation (Villumsen & Dav is|1986J : 



F(5) = 



1 + 



-1/2 



(19) 



The upper panel of Fig. [9] shows the distribution obtained com- 
puting the ratio between M raW|s and the real mass values M g ex- 
tracted from the simulation, for each galaxy in the FG subset (grey 
dots). The average and the la uncertainty region of the distribu- 
tion along the radial coordinate is also shown (thick solid line and 
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Table 2. Accuracy of mass estimation technique, evaluated by taking into account a mimicked observational error Auf in the line-of-sight velocity and the 
radial interval of aperture in observation. 
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[1.0,7.0] 


3693 


0.97 


1.46 


1.42 


0.67 





[1.0,3.5] 


1383 


0.96 


1.56 


1.49 


0.62 




[1.0,2.2] 


743 


0.93 


1.59 


1.48 


0.58 




[1.0,7.0] 


3746 


0.99 


1.48 


1.46 


0.67 


100 


[1.0,3.5] 


1421 


0.96 


1.58 


1.52 


0.61 




[1.0,2.2] 


752 


0.92 


1.62 


1.49 


0.57 




[1.0,7.0] 


3840 


1.02 


1.47 


1.50 


0.70 


200 


[1.0,3.5] 


1504 


1.01 


1.58 


1.59 


0.64 




[1.0,2.2] 


798 


0.97 


1.60 


1.57 


0.60 




[1.0,7.0] 


3761 


1.07 


1.51 


1.61 


0.71 


300 


[1.0,3.5] 


1488 


1.05 


1.62 


1.69 


0.65 




[1.0,2.2] 


802 


1.00 


1.67 


1.67 


0.60 



Note: /ijn and <tm are the moments of the distributions in the lower panel of Fig. [5] computed for different intervals of radial aperture. 



dashed lines, respectively). The radial coordinate is normalized to 
the virialization radius, rG,v,g = TG- g /r v . The distribution is sig- 
nificantly lower than 1 in the cluster core and approach unity in the 
non-equilibrium region. This is due both to projection effects and 
to the fact that the infall velocity approach is mostly reliable in the 
cluster outskirts, as already pointed out by CMM08. To correct the 
underestimation, we approximated the average distribution in the 
upper panel of Fig. [9] with a linear fitting algorithm and used it as a 
corrective term to compute a 'guess' mass value, as follows: 



Mg, 



M{?G,,; g )'- 



(20) 



where log 10 M (m, 



0.025f G ,v; 9 



0.26. 



The final corrected estimate of the cluster mass profiles, 
Af corr;9 , was obtained by smoothing the distribution of Mg ;9 val- 
ues for each cluster, taking its running median along the radial co- 
ordinate. We computed the running median within a window de- 
fined to contain approximately one-third of the galaxies in each 
cluster. The lower panel of Fig. [9] shows the distribution obtained 
computing the ratio between Af CO rr; 9 and the real mass values M g 
extracted from the simulation, for each galaxy in the FG subset 
(grey dots). Also in this case, the average and the la uncertainty 
region of the distribution along the radial coordinate is also shown 
(thick solid line and dashed lines, respectively). The combination of 
correction and running median smoothing yields very good results. 
The overall distribution is very close to unity from 1 to 7 virializa- 
tion radii. The variance of the distribution decreases when moving 
outwards, confirming once more the reliability of the infall velocity 
approach in the cluster non-equilibrium region. 

The overall accuracy of our technique was evaluated by com- 
puting the statistical moments of the distribution of M COII - g /M g 
values. We denote with (1m and om the logarithmical mean and 
the logarithmical variance computed among the set of values 
M corr /M: 



Hm = exp 



-x 



n ga i 



3 = 1 



In M cor 



InM, 



CM = 



exp 



1 / \ - 

— V} ( In M corr;9 - In M g - fx M J 
gai g=l V > 



(21) 



(22) 



ing into account the effects of the observational error and of the 
radial aperture in cluster observation. In particular, 

(i) we mimicked the observational error associated to the line- 
of-sight velocity measurement by substituting ui OS ; 9 with a per- 
turbed line-of-sight velocity ui os]g : 



^los; 9 — ^los;o ~\~ Au\ os 



id(0, Au^ 



(23) 



These moments were computed under different conditions, by tak- 



Here, rand(/i, it) is a function which generates random values hav- 
ing a Gaussian distribution with mean /i and variance a. We set 
Am!* os alternatively equal to 0, 100, 200 and 300 km s" 1 ; 

(ii) we defined different intervals of rG, v;9 , namely [1,7], 
[1,3.5], and [1,2.2], respectively, corresponding to the extent of 
uniform data coverage in our data catalogue, to the region encom- 
passed by turnaround radius (CMM08; CMM10), and to the tran- 
sition region between the cluster core and the non-equilibrium re- 
gion (where the infall approach is expected to be less reliable; see 
CMM08). The cluster core was rejected in all cases to avoid the 
contamination from virialization-related phenomena, which under- 
mine the reliability of infall approach. 

The values of hm and a M obtained with different choices of 
Aitfoj. and of the aperture interval are listed in Table [5] The ac- 
curacy is very good in all cases, confirming the reliability of the 
infall velocity approach applied to the FG subset. The variance o"m 
is only slightly affected by the observational error on the line-of- 
sight velocity and by the radial coverage in the galaxy distribution, 
confirming the robustness of the estimation. 

Our technique is also suitable to compute the virialization 
mass and turnaround mass of clusters. A/ corr , v and M corr ,t were 
obtained by interpolating the distribution of M corr;9 values at 
rG,v = 1 and at the radius where the overdensity equals the 
turnaround value S t = 15 (CMM08; CMM10). The interpolation 
was done with a linear fitting algorithm. In fact, this procedure is 
ineffective when the running median does not cover the radial ex- 
tent from the virialization radius to the turnaround radius. We de- 
note with /faii.v and /f a ii,t> respectively, the fraction of clusters for 
which the estimation of M CO rr,v and M corr ,t is impossible. Focus- 
ing on the remaining cases, we denote with /agree, v and / ag ree,t, 
respectively, the fraction of clusters for which the estimated values 
are consistent with the actual values A/ v and M t within the uncer- 
tainty. The values of / ta ii, v , /a grC c,v, /feil.t, and /a g rcc,t obtained 
for different choices of Aitj^ are listed in Table [5] The percentage 
of failure is slightly increased when the observational error in the 
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Table 3. Accuracy of the virialization and turnaround mass estimate, eval- 
uated by taking into account a mimicked observational error Am? in the 
line-of-sight velocity and the radial interval of aperture in observation. 



Au* oa (kmh- 1 ) 


jfail,v 


J 'agree, v 


/fail.t 


/agree, t 





0.10 


0.90 


0.06 


0.72 


100 


0.15 


0.93 


0.04 


0.69 


200 


0.19 


0.91 


0.03 


0.59 


300 


0.19 


0.91 


0.04 


0.54 



Note: /fail v an d /fail t are me fractions of clusters for which the 
estimation is impossible (computed among all clusters), /agree, v and 
/agree, t are the fractions of clusters for which the estimated value is in 
agreement with the actual value within the uncertainty (computed among 
the clusters for which the estimation is possible). 



line-of-sight velocity is introduced, as an effect of the higher noise 
in the galaxy velocity distribution. Nevertheless, the estimation of 
both A/ v and A/ t is significantly accurate even for relatively high 
values of Au* os . 

4.1 Single-cluster estimations 

Some characteristics of our mass estimation technique can be bet- 
ter appreciated by comparing the mass profiles estimated for single 
clusters separately. In this section, we will discuss in detail the re- 
sults obtained for nine simulated clusters from our data catalogue 
and for nine observed clusters from the CIRS Catalogue by RD. 

Fig. [10] shows the mass profiles of the nine simulated clus- 
ters. The grey solid lines represent the actual profiles M g , while 
the black solid lines and dashed lines represent the estimated profile 
M COII - g , together with its uncertainty band. We adopted the 2<r un- 
certainty on the running median as the uncertainty on this estimate, 
since it was observed to mimic well the dispersion in the actual 
distribution of Mg-, 3 values. The distributions of values M r aw; S are 
also shown with grey dots. This subset of nine clusters gives a good 
representation of the common results obtained when estimating the 
mass of single clusters from the whole catalogue. In most cases, 
our technique is able to produce a mass profile extended from the 
virialization core up to the extreme outskirts of clusters (clusters 1, 
2, 4, 6, 7 and 9), even when only a few FGs are identified (clusters 
1 and 7). If the FGs are not evenly distributed, the estimated profile 
turns out to be radially limited (clusters 3 and 8) but the estimate is 
generally correct within the covered interval. The definition of the 
FR, which is the same for all clusters (as described in Section [3T) 
does not mask the individual features of different clusters, which 
are correctly reconstructed through the features of the FG subset. 
In fact, our technique is able to discern the presence of remarkable 
changes of slope in the mass profiles (clusters 1, 2, 3, 6 and 9), 
which cannot be predicted using the results of CMM08. In some 
cases, the procedure of correction and smoothing produces an un- 
desired plateau in the outermost regions (clusters 9). Only in few 
cases, the poor quality of data prevents our technique from properly 
reconstructing the mass profile (cluster 8). 

The values of the virialization mass and the turnaround mass 
estimated for the nine simulated clusters in Fig.[K)]are listed in Ta- 
ble|4] along with the corresponding values extracted from the simu- 
lations. In most cases, the agreement is fairly good, even when the 
overall profile is poorly reconstructed, thanks to the linear fitting 
algorithm used to compute A/ cor r,v and M COTT>t . 

Fig. [10] shows the mass profiles of the nine CIRS clusters. 
These clusters were chosen among the ones with the highest ra- 



dius of complete radial coverage in the survey. As before, the black 
solid lines and dashed lines represent the estimated profile A/ CO rr; 9 , 
together with its uncertainty band, while the distributions of values 
M r aw; 9 are shown with grey dots. The virialization radii r v of these 
clusters, needed to normalized the galaxy radial distribution, were 
computed iteratively as the radii where the overdensity estimated 
by our technique equals the virialization value S v ~ 101/Slo — 1 
I Bryan & Norman|l998] CMMOS^The CIRS clusters are gener- 
ally more populated than those extracted from our simulation, al- 
lowing a lower uncertainty in the mass profile estimation. A fairly 
good agreement with the estimates of RD is generally observed at 
all radii. Compared to the results of RD, our estimate is generally 
less accurate in the cluster core and better accurate in the outskirts. 
In most cases, our technique is able to reconstruct the mass profiles 
well beyond the turnaround radius. 

The values of virialization mass and turnaround mass esti- 
mated for the nine CIRS clusters in Fig. [TTJ are listed in Table [5] 
along with the corresponding values estimated by RD. Only in 
this case, we adopted as virialization mass the value M cor r,200, 
which corresponds to the commonly used overdensity value S200 = 
200/fio — 1, to allow better comparison with the values A/200 es- 
timated by RD. We point out that such estimation of Af CO rr,200 is 
meant only to provide a reference value and must not be regarded 
as reliable, because r2oo is smaller r v and falls outside the effec- 
tiveness interval of our tecnique. Nevertheless, the values of M200 
and A/ C orr,200 are in quite good agreement within the uncertainties. 
A fairly good agreement is observed also between A/ t and M corrjt . 
In some cases (clusters A971, A1436, A1885) our estimates of the 
turnaround mass are significantly higher than those of RD. This is 
not unexpected, since the profile of these clusters as estimated by 
RD are truncated before actually reaching r t , resulting in a possible 
systematic underestimation of A/ t . 



5 DISCUSSION AND CONCLUSION 

We discussed a new approach to estimate mass profiles and masses 
of galaxy clusters, based on the relation existing between the cluster 
overdensity and the infall velocity of member galaxies ( |Regos &] 
|Geller|1989} . Our technique is simple, it only needs to know sky- 
plane positions and the line-of-sight velocities of the galaxies, and 
it can be applied to observed clusters. 

Our analysis was performed and checked on a simulated cat- 
alogue of clusters (Borgani et al. 2004, Biv iano et al.|2 006) in the 
theoretical framework of the SCM which is widely accepted in lit- 
erature (CMM10). We found the existence of an FR in the redshift 
space and we demonstrated that the galaxies belonging to this re- 
gion, i.e. the FGs, are able to effectively identify the mass profiles 
of clusters and to measure the corresponding total masses. 

Our technique consists in: 

(i) Identification of the FG subset via the identification of the 
FR in the redshift space; 

(ii) Detection of the mass profiles using the redshift-space coor- 
dinates of FGs. 

1 An alternative approach for computing r v is provided by the SCM. 
According to the results CMM08 and CMM10, the ratio between the 
turnaround radius and the virialization radius is approximately constant: 
r t /r v ~ 3.5. Therefore, r v can be obtained from r t , defined as the ra- 
dius were the estimated overdensity equals the turnaround value <5t ~ 15. 
This approach is suitable for estimating the values of r v and M v in poorly 
populated clusters, when no FGs are identified within the virialization core. 
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Figure 10. Mass profiles for nine simulated clusters from our data catalogue. In each panel, we represented the actual cluster mass profile M s (grey solid 
lines), the distribution of M raw;g values (grey dots) and the extimated mass profile M CO rr;g (black solid lines) together with its uncertainty band (black dashed 
lines; see Section[4}. The turnaround radii rt are also shown (vertical dotted lines). 



Table 4. Estimated virialization masses and turnaround masses for the nine simulated clusters in Fig. Qo] compared with the corresponding values extracted 
from the simulation. 



cluster M V (10 14 M Q ) M corr , v (1O 14 M ) Aft (1O 14 M ) M cor r,t (1O 14 M ) 



1 


1.37 


i H7 +2.69 
L -°'-0.8S 


2.48 


2.86t^ 7 


2 


1.75 


i-8oi :lS 


3.87 


s.ostiS 


3 


2.07 




8.13 


9-24tr 4 8 9 


4 


3.36 


1 98+ 189 


7.04 


6 44+ 1 74 


5 


3.64 


, 99 +1.17 
°- z -0.59 


6.40 


7 19+ 244 


6 


4.08 


o co+0.89 
o.oo_ 114 


8.25 


_ --+0.98 
'•' '-1.87 


7 


1.16 


1 QQ + 1-58 

i-l«_0.41 


5.26 


5 n 7 +5.46 
°- 9 ' -0.75 


8 


1.17 


Z - ZO -0.89 


3.45 




9 


4.11 


A 40+2.33 
4-*0_ 1 _ 47 


7.12 


D -^°-0.71 



We showed that it is possible to estimate the cluster mass profiles 
from 1 up to 7 virialization radii, within a tipical uncertainty fac- 
tor of 1.5, for more than 90 per cent of clusters. The technique is 
reliable even with few identified Fair galaxies. 

Our technique was tested to be accurate for clusters in the 
mass range provided by our simulation, corresponding to virializa- 
tion masses between 8.0 x 10 13 Ii^Mq and 1.3 x 10 15 h~ 1 M (s . 
In the present paper, we chose for the cosmological matter density 
parameter the value Qm, = 0.3, which is also the value adopted 



by the simulation we used. Our results can be adapted to different 
values Om.o. ranging from 0.2 to 0.4. 

We applied our technique in a subset of clusters taken from 
the CIRS Catalogue (RD). The mass profiles and the masses we 
obtained are in fairly good agreement with previous literature val- 
ues. 

In the future, we aim to apply the present technique to clusters 
of large observational catalogues. 
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Figure 11. Mass profiles for nine observed clusters from the CIRS Catalogue (RD). In each panel, we represented the distribution of Af raw;g values (grey 
dots) and the extimated mass profile M corr;9 (black solid lines) together with its uncertainty band (black dashed lines, see section]?). The virialization radii 
T2oo and the turnaround radii rt are also shown (vertical dotted lines, leftmost and rightmost in each panel, respectively). 

Table 5. Estimated virialization masses and turnaround masses for the nine CIRS clusters in Fig. |l 1| compared with the corresponding values obtained by RD. 



cluster M 2OO (1O 14 M ) M corr , 2 oo (1O 14 M ) M t (1O 14 M ) M corr , t (1O 14 M ) 



A671 


3.23 ± 1.02 


1 ™+ - 82 

l.OU_ 29 


4.40 ± 1.50 




A971 


4.46 ± 0.89 


2 Qfi+ 2 15 
Z - MD -1.81 


4.95 ± 1.05 


io.4±?;H 


A1066 


4.68 ±0.51 


1 98+ ' 49 
1 ' M °-0.91 


7.93 ± 1.06 


6 25+ 102 


A1314 


1.72 ±0.44 


2 42+ 104 
-0.39 


2.28 ±0.69 


*- MD -0.78 


A1436 


0.86 ±0.22 


, 07 +1.45 
-2.91 


1.83 ±0.67 


c 7 9 +1.24 


A1767 


6.03 ±2.36 


2 4fi+ 1;33 


11.3 ±4.88 


fi 7fi +1.56 
D -' D -0.65 


A1773 


1.98 ± 1.90 


, 79 +0.73 
' -0.88 


3.30 ± 3.18 


4 98+1-01 


A1885 


4.50 ± 1.24 


2 59+ 4 33 
Z - Ja -0.81 


4.46 ± 1.23 


6.82i 2 ;^ 


A2064 


1.65 ±0.63 


0.23tlit 


4.27 ± 1.48 


c -JQ+0.68 
O-jy_0.87 
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